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Jupiter's composition: sign of a (relatively) late formation in a 
chemically evolved protosolar disk 
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ABSTRACT 

It has been proposed that the enrichment in noble gases found by Galileo in Jupiter's atmo- 
sphere can be explained by their delivery inside cold planetesimals. We propose instead that 
this is a sign that the planet formed in a chemically evolved disk and that noble gases were 
acquired mostly in gaseous form during the planet's envelope capture phase. 

We show that the combined settling of grains to the disk midplane in the cold outer lay- 
ers, the condensation of noble gases onto these grains at temperatures below 20-30K, and the 
evaporation from high disk altitudes effectively lead to a progressive, moderate enrichment of 
the disk. The fact that noble gases are vaporized from the grains in the hot inner disk regions 
(e.g. Jupiter formation region) is not a concern because a negative temperature gradient pre- 
vents convection from carrying the species into the evaporating region. We show that the ~ 2 
times solar enrichment of Ar, Kr, Xe in Jupiter is hence naturally explained by a continuous 
growth of the planet governed by viscous diffusion in the protosolar disk in conjunction with 
an evaporation of the disk and its progressive enrichment on a million years timescale. 

Key words: Solar system: formation, planetary systems, planetary systems: formation, plan- 
etary systems:protoplanetary disks, accretion disks, planets: Jupiter, planets: Saturn 
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1 INTRODUCTION 

Ten years have elapsed since the Galileo probe entered Jupiter's 
atmosphere and allowed in situ measurem ents of its composition 
(e.g.. lYoung et alll996tlOwen et all! 999). One of the key results 
of the mission is that, when compared to hydrogen, most measur- 
able species appear to be enriched by a factor two to four compared 
to a solar composition. 

Most surprisingly, the list of enriched species includes not 
only "ices" (i.e. species that condense at temperature ~ 100 K 
in the protosolar disk like methane, ammonia), but also the noble 
gases Ar, Kr, Xe tMa haffv et aD boOO). which are enriched com- 
pared to the new solar abundances (e.g. Lodders 2003, for a re- 
vie w) by factors 2.2 ± 0.8, 2.0 ± 0.7, 2.0 ± 0.5, respectively. This 
led lOwenetail {l999) to propose that Jupiter formed from very 
low-temperature planetesimals, at a temperature ~ 30 K allowing a 
direct condensation of noble gase s onto a morphous ice. This expla- 
nation was challenged by Gautier et aL 1 2001) who proposed that 
noble gases would be incorporated into crystalline ice by clathra- 
tion at slightly higher temperatures, and that they would thus be 
delivered into Jupiter by icy planetesimals. The large amount of 
cages available to trap noble gases implies that the ratio of oxy- 
gen to other species in Jupiter should be highly non-solar, but the 
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precise factor depends on the condensation se quence iGautier et alj 
l200ltlHersant et all20ollAlibert et all2005l) . 

In this letter, we propose instead that the present abundances 
of Jupiter indicate that the planet (and other giant planets as well) 
formed in an chemically evolved protoplanetary disk, in which part 
of the light gases (hydrogen, helium and maybe neon) had been lost 
by photoevaporation. 

Except as oth erwise noted, we will use the disk evolution 
model described by H ueso & Guillotl i2"b05) with the following val- 
ues of the parameters: M c a = 1 Mq (mass in the molecular cloud 
core), Mq = 0.4 Mq (seed mass for the protosun), T c( j = 10 K 
(ambient interstellar temperature), o; c d = 10 _14 s _1 (assumed uni- 
form rotation rate of the cloud core), a — 0.01 (viscosity). 



2 SOLIDS AS CARRIERS OF THE NOBLE GASES? 

Explaining the efficient delivery of noble gases into Jupiter's at- 
mosphere by embedding them into solids (either through a direct 
condensation at ~ 30 K, or through clathration at slightly higher 
temperature) is difficult: As shown by Fig.Q the disk should al- 
ways remain warmer due to stellar irradiation up to 10 AU from 
the central star, i.e. where Jupiter and Saturn should have formed. 
A possibility could be that the disk self-shadows this region (e.g., 
iDullemond et alj 2001). but it remains to be proved that this effect 
occurs at the right time and distances, and that it is able to maintain 
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Figure 1. Temperature pressure diagram showing the region of condensa- 
tion of H2O, and the clathration lines for Xe, Kr and Ar into crystalline 
ice, assuming a solar composition. The coloured lines are evolution tracks 
in the protosolar disk at distances of 5, 10, 25 and 100 AU, respectively, 
for irradiated disks (plain) or without including stellar irradiation (dashed). 
The symbols corresponds to ages of 1 (squares), 2.5 (diamonds) and 5 Ma 
(triangles), respectively. The disk evolution models are calculated without 
evaporation. 

sufficiently low temperatures. Similarly, Jupiter could have formed 
further from the Sun and migrated inward 1 Alibert, Mousis & Benz 
2005), but it is yet unclear why Jupiter would have accreted its gas 
in the > 10 AU region, sufficiently late for the disk to have cooled, 
and then not in the 5 — 10 AU region when its migration there im- 
poses that a massive disk was still present. 

Furthermore, the clathration model requires a large number 
of available cages to trap a significant fraction of the argon. This 
implies that the ratio of H2O to all other clathrated species has 
to be large to provide as many cages as possible, and yet it is 
constrained by interior models of Jupiter that imply that the to- 
tal mass of heavy elements in Jupiter, including the core and in- 
cluding rocks, is smaller than 40 M e 1 Guill otkOoH) . Alibert et al. 
explain the observed abundances of noble gases by the accretion 
of ~ 20 Me of clathrated ices, assuming a 100% clathration ef- 
ficiency. An efficiency of 50% or less would require accreting a 
larger amount of ices, in violation of the interior models constraints. 

In summary, although it appears very likely that chemical 
species including noble gases were at some point trapped into solid 
planetesimals, explaining Jupiter's composition through that mech- 
anism only seems to require an unlikely combination of factors. 



3 GIANT PLANET FORMATION IN A CHEMICALLY 
EVOLVED DISK 

We now assess the possibility that the giant planets were formed 
in a chemically evolved protoplanetary disk and that noble gases 
were delivered to the planets mostly in gaseous form. In our model, 
the noble gases are trapped into the solids at low temperatures in 
the outer disk, but most of them are released in gaseous form in 



the giant planet formation region. We postulate that this trapping is 
almost complete at 30 AU and beyond because of the low (~ 20 K) 
temperatures (see lNotesco & Bar-Nunl2005l and fig.0. 

For simplicity, we consider that our protoplanetary disk is 
made of 2 components: a dominant one that includes hydrogen, he- 
lium and neon (species that remain gaseous at all temperatures), of 
abundance xu ~ 1; and a minor one (e.g. argon, krypton or xenon) 
of abundance xz = x <C 1 that is carried with the grains to the 
midplane of the disk and radially inward. 

3.1 Enrichment by an inward radial migration of grains? 

The inward radial migration of grains of ~cm size has been postu- 
lated to be a key ingredient in the problem of the formation of plan- 
etesimals (e.g. ICuzzi & Zahnld 12004). Since these grains should 
have captured noble gases in the cold outer regions of the disk, 
this could potentially also lead to an enrichment of the inner disk in 
noble gases. We find however that this scenario is unlikely to yield 
noble gases abundances that are compatible with the measured val- 
ues in Jupiter. This is due to two reasons: 

(i) Given that we expect a fast mixing of the inner regions, the 
local enrichment of a given species should be proportional to the 
disk mass inside its vaporization radius. Now, this factor is very de- 
pendent on the vaporization temperature: Compared to a species 
with a 20 K vaporization temperature, species that would be re- 
leased at 30, 40 or 50 K should be more enriched by factors ~3, 
5 and 10, respectively. Given the relatively uniform enrichment 
of Jupiter's atmosphere in Ar, Kr, Xe, and given that Ar is much 
more volatile than Xe, it is difficult to imagine that this scenario 
can work. 

(ii) An efficient transport of noble gases into the disk's inner 
region also implies a rapid loss by accretion onto the star, and a 
global depletion of these elements. Solving the problem then re- 
quires a careful balance between inward transport and loss of the 
depleted outer disk (which otherwise will be accreted onto the giant 
planets). 

Preliminary mass balance calculations indicate that the inward 
transport of noble gases had to be very limited, and is unlikely to 
have had an important role: we choose to neglect it in what follows. 
We will assume that it may have occurred early on, but led to the 
loss of only a minor fraction of the noble gases in the disk. 

3.2 Evaporation from the disk 

3.2.1 Main characteristics 

It is relatively natural to invoke evaporation of disk material in order 
to explain the loss of light hydrogen and helium while retaining 
heavier elements. The fact that disks evaporate has been observed 
in the relatively extreme situations of proplyds under the intense 
irradiation of O-B stars, but is also commonly invoked to account 
for the fast disappearance of disks around young stars. 

The temperature in the disk atmosphere controls the evapora- 
tion, and defines the critical radi us, when t he sound speed is equal 
to the orbital velocity (e.g. Shu et al. 1993): 



where G and k are the gravitational and Boltzmann constants, M* 
is the mass of the central star, and fj, the mean molecular mass. For 
a solar- type star and atomic hydrogen, R 3 is of order of 10, 100 and 
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Table 1. Average characteristics of the evaporation flow estimated at the disk's mid-life 
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300 AU, for atmospheric temperatures of 10000, 1000 and 350 K 
respectively. 

It is generally thought that only disks that are supercritical 
(i.e. r > R a somewhere i n the disk) evaporate l Shu et al. Il993t 
iHollenbach et alfeOOCtlciarke et alfcOOll) . but lAdams etalJCOOfl) 
show that an efficient evaporation is also possible in subcritical 
disks. In all cases, this requires that a disk atmosphere is super- 
heated significantly compared to the midplane. High temperatures 
(~ 10 4 K) in the inner disk ( ~ 10 AU) can be due to extreme UV 
radiation from the central star I S hu et aljl993l) : Moderate tempera- 
tures (~ 100 — 600 K) in the outer disk can be maintained by the far 
UV ambient radiation from a stellar cluster lAdams et all2004 or 
by disequibrium heating of the ga s in the presence of som e dust par- 
ticles ( Jonk heid et all2004lKamp & Dullemonj2004) . (Far more 
violent evaporation effects can occur near massive stars, but for 
simplicity, we choose not to discuss this possibility). 

Table Q depicts the global characteristics of the evaporation 
flow. A crucial point is that the evaporation is hydrodynamical, i.e. 
it occurs at levels where the mean free path A csc is significantly 
smaller than the pressure scale height Hp. Because of the very low 
gravity, the evaporation also occurs without any separation of the 
chemical elements, as shown by the high value of the critical mass 
(Huntenet all 1987b: 



rit/mH = 1 + 



fcTE c 
bgXHm^ ' 



(2) 



where & is a diffusion coefficient (we use b ~ 10 cm 
the gravity. 

In all cases, the levels at which evaporation takes place cor- 
respond to relatively high altitudes in the disk, much higher than 
the estimated thickness of the dust su bdisk I Dubrulle^UilJ^95j; 
[ Dullemond & Dominildl2005t iTanaka et all2005t iThroop & Ballvl 
2005). This implies that species that condense onto grains and are 
transported towards the midplane will generally not be lost from 
the disk. Furthermore, because a negat ive vertical temperature gra- 
dient prevails (e.g. iMalbet & Bertoutlfl99lt IChiang & Gol dreich 
Il997t) . even species that are captured by grains at low temperatures 
but later vaporize when subject to higher temperatures will not be 
transported by convection back up to the levels where they would 
hydrodynamically escape. This implies that the disk can become 
progressively enriched in species that condense onto grains at low 
temperatures down to ~ 15 — 30 K. 



3.2.2 The model 

Let us consider the evolution of a viscous accretion disk of sur- 
face density E(r, t), subject to an evaporation ra te S eBC ■ Its evolu- 
tion is governed by the following equation (e.g. lciarke et all200lt 
iHueso & GuillofcOOSh : 



dt 



3d_ 

r dr 



,1/2 



l/2\ 



(3) 



and the escape rate is limited to a layer of temperature T csc , mean 
molecular mass m CS c, corresponding sound speed c CS c, and number 
density n csc : 

Ecsc W'escCesc^'esc • (4) 

Let us now consider that a second minor component of surface 
density cE is present (with c <C 1) and obeys the same diffusion 
equation as for E, with the exception that for the reasons discussed 
previously, its escape rate is very small. (Obviously, this is a sim- 
plification because one should include backreactions between the 2 
species). Within this framework, the evolution of the concentration 
c, or equivalently of the enrichment £ = c/cq (where cq is the 
protosolar concentration of that element) can be shown to obey the 
following equations: 



d£ , 



£ =11- 

t-'CSC I J- 



, dln(i/£) 
<91nr 



3\ ld£ 
+ 2 I r dr + 



d 2 £ 

dr 2 



En 



(5) 



(6) 



Here, we assume that the enrichment can be at most £ max , 
after which both the major and the minor species escape. This 
maximal enrichment can be estimated as follows: (i) If the evap- 
oration takes place inside the noble gases vaporization radius, the 
hydrogen-helium atmosphere of the disk is progressively eroded 
until the midplane conditions are reached. This implies values of 
£max of the same order as the midplane dust enrichment estimated 
at the vaporization radius. Depending on properties of the d ust and 
turbu l ence, it could be eit her large (<; 10) (Dullemond & Dominik 
l2005t ITanaka et al]|2005h or small (< 2) <Dubrulle et alJll995h . 
(ii) For an evaporation of the outer disk, the condensation of no- 
ble gases onto grains implies that the scale height of the con- 
densing species h is maintained at a smaller value than that of 
the gas, H. Because the vertical density profiles are gaussian, 
the maximum enrichment can be shown to be of order £ max ~ 
H/hevfc(z eB c/Hp)/evk(z eB c/h), where z csc is the altitude at 
which escape takes place. The fact that z csc > H in most cases 
garanties that f max is large (i.e. 3> 10), even for values of H/h 
close to unity. Globally, our model rests on the likely assumption 
that <f max is large, and we use a fiducial value £ ma x = 10. 

Note that equation j5J contains both a diffusion term and an 
advection term. The latter accounts for the fact that both near the 
inner and outer boundaries, matter is essentially advected as it is 
either accreted by the star or lost from the system, respectively. Ba- 
sically, the problem to be studied is one in which the disk becomes 
progressively enriched in heavy elements as it loses hydrogen and 
helium, but at a rate that is limited by the diffusion of these heavy 
elements in the disk, with a timescale ~ R 2 /v. 
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We further assume that the accretion of the envelopes of 
Jupiter and Saturn is limited by the viscous diffusion in the disk 
so that the planets' growth is Csee lLecar & Sasselovl2 003'): 



M p j a net (£) = e p laiict A/disk (^planet, t) , 



(7) 



where MdiskfVpianet, ' s m e mass flux in the disk that crosses the 
annulus of radius r p ianet at time t. The parameter e p i an et is the ra- 
tio of the flux accreted by the planet to the total flux that visously 
diffuses inward. It is mostly unknown, and certainly very depen- 
dent on the mass of the planet and its cooling. Based on numerical 
simulations for Jupiter, we choose e p ianet ~ 0.3 (F. Masset, pers. 
communication). 



3.2.3 Inner disk evaporation 

We first follow the approach of Clar ke et"all 1200 ll) by assuming 
that the young T-Tauri Sun emits extreme UV photons at a rate $, 
and expecting that <E> « 10 41 s _1 . In that case, given that the pho- 
tons heat the upper layers of the disk to ~ 10, 000K, the evapora- 
tion m ostly takes place at a critical radius R g ~ 10 AU IShu et ail 
119931) . The escape flux is defined by eq. with T csc « 10 4 K and 
by 



n csc = no$ 
n csc = 



1/2 rv 
41 "-gli 



-5/2 



for r ^ R g 
otherwise 



(8) 



where n = 5.7 x 10 4 cm" 3 , $41 = $/10 41 s _1 and R gli = 
J? 9 /10 14 cm. The fact that viscous diffusion is much fa ster close to 
the Su n leads to an inside-out removal of disk material l Clar ke et"all 
1200 lh . 

Figure [2] shows the resulting evolution of the disk mass and 
disk enrichment, the growth of Jupiter and Saturn as derived from 
the disk properties and eq. Q and the consequent enrichment of 
their envelopes, assuming that they are fully mixed by convection. 
The effect of evaporation becomes apparent only when mass of the 
disk decreases to a point when the mass flux in the disk becomes 
comparable to the evaporation mass flux. Because the evaporation 
proceeds from inside out, inner regions (< 10 AU) are depleted 
rapidly and the corresponding local enrichment (2nd panel in fig.|2j 
can become very large. When this occurs, the mass flux accross 
these radii has become tiny, so that Jupiter's enrichment is always 
found to be moderate and close to the observed values. 

We find that Saturn should be about 50% more enriched ac- 
cording to this scenario. However, this crucially depends on the 
assumption that Saturn and Jupiter stole the same constant part of 
the disk mass flux to build their envelopes (i.e. ej up == es a t = 0.3). 
If Saturn has been less efficient at that, then it may be as enriched 
in noble gases as Jupiter. 



3.2.4 Outer disk evaporation 

There are two problems with the EUV evaporation scenario: (i) The 
timescale for the loss of the circumstellar disk is uncomfortably 
long compared to observations, even in the extreme case of $ = 
10 43 s _1 ; (ii) It requires a high, constant value of EUV photons 
production fro m an u nidentified mechanism not powered by the 
accretion flow (Matsuvama et al. 2003). 

Adam s et alJ 1200 4) hence proposed a scenario based on a 
moderate heating of the disks atmosphere by ambient FUV radi- 
ation and a subcritical evaporation. We use a simplified version 
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Figure 2. Time evolution of the mass of the protosolar disk Mjiak. of its 
local (5 AU) enrichment fdisk. °f tne masses Mpi ane t and enrichments 
of Jupiter (plain) and Saturn (dashed). The dotted portions of the 



planet 



curves correspond to the assumed (slow) growth phase of the planets' cores. 
The calculations are made for a EUV evaporation flux of 10 41 (blue) and 
10 43 s — 1 (red), respectively. The dashed area in the bottom panel indicates 
the enrichments in Ar, Kr and Xe inferred for Jupiter from the Galileo probe 
data. 



of their evaporation rates (see their appendix) to obtain an evap- 
oration that is controlled by the temperature of the escaping disk 
atmosphere: 



EUV I ih. \ 1/2 p -(R S /r-l)/2 



forr ^ R g 
otherwise, 



(9) 



where ctpuv ~ 10 -21 cm 2 is the FUV dust cross section. These 
escape rates correspond to the evaporation from the disk surface. 
The mass lost in the radial direction is to be considered and is even 
dominant for subcritical disks, and we hence remove the mass to 
the last radial layer i?d of the disk: 

Mrad = rn-csc 

_1 / R \ ' 



for R d > Rg 

-(«s/«d-l)/2 



(10) 



otherwise 
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Figure 3. Same as fig.lSI but in the case of an evaporation of the outer disk 
for three temperatures of the escaping disk atmosphere: 100, 200 and 600 K. 



Figure [3] shows the evolutions of the disk and planets result- 
ing from the combination of early supercritical and later subcritical 
evaporation, as the disk progressively shrinks. As radial evapora- 
tion does not lead to an increase in £ , the enrichment is always 
self-limited, and mostly controlled by the value of T csc (i.e. by the 
global ratio of the vertical to radial evaporation). Once again, we 
find that for a quite wide range of values of the different parame- 
ters, the resulting enrichment for Jupiter agrees with the measured 
enrichment in noble gases. In this scenario, because the enrichment 
in the disk is very progressive, it is expected that Jupiter and Saturn 
should have about the same enrichment in noble gases. 

Even though the details of the photoevaporation have been 
greatly simplified, we stress that this model does provide an ex- 
planation for the Ma timescales for the removal of circumstellar 
disks, with viscous diffusion playing a key role in spreading the 
disk outward to regions where it can evaporate more efficiently. 



3.2.5 Sensitivity to input parameters 

We have shown that different scenarii of the evaporation of disks 
lead to similar results in terms of final enrichments of the giant 



planets. We briefly describe how the most important parameters 
affect the results: 

Viscosity (a): It affects the global evolution timescale (pro- 
portional to a in the outer-disk evaporation scenario), but very 
marginally the final values of the enrichment. 

Initial angular momentum (u) c d)'. The maximum disk mass is 
directly related to the amount of angular momentum of the molec- 
ular cloud core from which the disk is born. We find that except for 
low values of this parameter, the final £ p i an ct remains quite similar. 

Maximal enrichment (£ ma x): Our model rests on the assump- 
tion that the local enrichment in the disk can be relatively large; the 
conclusions of the article remain qualitatively valid and quantita- 
tively similar as long as £ m&K 5. For smaller values we expect 
the final planet enrichment to be smaller than calculated here and 
incompatible with the observations. 

Planet growth factor (e p i anct ): This parameter governs the 
growth of our model planets. Because the enrichment in the disk 
rapidly increases, larger values of this factor yields larger £ p i ane t- 
However, the dependance is only moderate because of the rapid de- 
cline in the disk (and planet) accretion rate. We tested that changes 
are smaller than the error bars on the abundances of noble gases for 
values of e between 0.1 and 1. 



4 CONCLUSIONS 

We have shown that the observed enrichment of Jupiter in noble 
gases can be explained by a progressive capture of the planet's en- 
velope in sync with the evaporation of the protosolar disk. We in- 
fer that the envelope capture phase started relatively late, at about 
half of the disk's lifetime, probably because of a slow or delayed 
growth of the protoplanetary core. In our scenario, the final masses 
of Jupiter and Saturn result naturally from a competition between 
limited viscous accretion and disk evaporation. We note that a wel- 
come consequence of the late start scenario is a significant sup- 
pression of the inward migration when compared to models in 
which giant planets form early, in locally massive disks (such as the 
minimum mass solar nebula or more). Provided disk atmospheres 
can be heated to temperatures of 100 K or more, the disk model 
that we propose further explains the disappearance of disks in Ma 
timescales, as observed. 

We predict that the enrichments in Ar, Kr, Xe in Saturn should 
be similar to that in Jupiter. It should be larger in Uranus and Nep- 
tune, due to a late capture of more enriched gas from the evaporat- 
ing disk, but also to a pollution of unknown magnitude from their 
ice cores, significantly more massive than their gaseous envelopes. 
No isotopic fractionation is expected to arise from the disk evapo- 
ration process. However, clues on the ambiant mean temperature of 
the protosolar molecular cloud core may be obtained from the pri- 
mordial Ne/Ar, Ne/Kr, and Ne/Xe ratios: if this mean temperature 
was high enough (e.g. 20 ~ 30 K) not to allow neon to condense 
onto grains, then neon should have escaped progressively from the 
disk with hydrogen and helium, yielding small values of these ra- 
tios. In all cases these ratios should be equal to or smaller than 
unity. 

These results show the importance of noble gases for tracing 
back events that occurred in the early Solar System and stress the 
need for an accurate determination of the compositions of all giant 
planets. 
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